High-Energy Limit of Massless Dirac Fermions in Multilayer Graphene using 
Magneto-Optical Transmission Spectroscopy 
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We have investigated the absorption spectrum of multi-layer graphene in high magnetic fields. The 
low-energy part of the spectrum of electrons in graphene is well described by the relativistic Dirac 
equation with a linear dispersion relation. However, at higher energies (> 500 meV) a deviation 
from the ideal behavior of Dirac particles is observed. At an energy of 1.25 eV, the deviation from 
linearity is ~ 40 meV. This result is in good agreement with the theoretical model, which includes 
trigonal warping of the Fermi surface and higher order band corrections. Polarization-resolved 
measurements show no observable electron-hole asymmetry. 

PACS numbers: 78.67.-n, 73.21.-b 
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Graphene, a single sheet of graphite, is a two- 
dimensional system which exhibits unique electronic 
properties mostly related to its peculiar band struc- 
ture P, 0, 0, 13] • The remarkable physics exhibited by 
graphene has its origin in the conduction and the valence 
bands which meet at the two inequivalent (K and K') 
corners of the Brillouin zone. The electrons in the vicin- 
ity of the Fermi energy do not obey Schrodinger's equa- 
tion, but should instead be described using the quantum- 
electrodynamic Dirac equation for relativistic fermions 
with zero rest mass. The electrons have a linear dis- 
persion relation whose slope defines a Fermi velocity vf- 
In a relativistic analogy, these electrons behave as mass- 
less Dirac fermions moving at an effective speed of light 
Vf- This system is of great interest from a fundamental 
physics point of view and it has even been suggested that 
graphene can be used for bench top quantum electrody- 
namics experiments for example to test the Klein 
paradox [5j. However, in graphene, considering the car- 
riers as massless fermions remains an approximation and 
it is both important and interesting to verify the limits 
of this approximation. 

Graphene has been extensively investigated using op- 
tical measurements such as Raman scattering [y, |7|, [8|, l2| , 
far-infrared absorption (FIR) 0,[lH, as well as magneto- 
photoconductivity [13]. Landau level (LL) spectroscopy 
is a direct and precise tool to test the linear dispersion 
relation in the close vicinity of the K and K' points of 
the Brillouin zone. In the presence of a magnetic field 
B, perfect linearity leads to the observed V Bn spacing 
for the LLs indexed by the integer n. In the low-energy 
range of the Dirac cone, the linearity of the dispersion re- 
lation is well preserved [13, [HI EH ■ However, graphene is 
a solid-state system composed of carbon atoms arranged 



in a honeycomb lattice and the linear dispersion in these 
specific high symmetry points is only a part of a com- 
plicated band structure. This implies that the analogy 
with neutrinos, massless Dirac particles, cannot hold ev- 
erywhere, and we expect a deviation from a linear dis- 
persion for high energies of the Dirac cone. Although 
it is well established that at low energies, electrons in 
graphene can be treated as massless Dirac particles, it is 
crucial to determine the limits of this approach. 

In this Letter, we probe the limits of the massless Dirac 
fermion approximation in graphene by extending the pre- 
vious studies [niEil to higher magnetic fields, and, most 
importantly, to higher energies. Using magneto-optical 
transmission spectroscopy, we present a full LL spec- 
troscopy in magnetic fields up to 32 T, from the far in- 
frared to the visible range of energy. Transmission mea- 
surements performed in the near visible provide an access 
to the high-energy range (< 1.25 eV) of the Dirac cone. 
A significant deviation from the linear dispersion of ideal 
Dirac fermions is observed. The experimental data are 
compared to a theoretical model which includes higher- 
order band terms and a good agreement is obtained. In 
addition, the asymmetry between electrons and holes has 
been probed using polarization-resolved transmission ex- 
periments. 

We have investigated samples containing a high num- 
ber of graphene layers (between 70 and 100) grown in 
vacuum by the thermal decomposition method, on a (4H) 
SiC 14|, ll5( substrate. Both experiment and theory con- 
firm that the layers are electronically decoupled so that 
the system can be considered as a multi-layer graphene 
sample 0, E(l El, EH] • In particular, these samples show 
Raman spectra with the characteristic signature of sin- 
gle layer graphene It is likely that this peculiar- 
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FIG. 1: (color online) (a) Representative differential trans- 
mission spectra (the B = T spectra has been subtracted) 
measured at the given magnetic fields, (b) Positions of the 
absorption lines as a function of the square root of the mag- 
netic field. Stars represent data obtained in near visible range, 
circles denotes the data measured by FTS. Dashed lines are 
calculated energy of the transitions between LLs assuming a 
linear dispersion. On the right hand side the observed tran- 
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and determined a Fermi velocity vp = 1.1 x 10 6 m/s. This 
is somewhat larger than the value Dp — 1.02 x 10 6 m/s 
obtained here from the slope in the low-energy region 
of Fig. Hlb). Moreover, in Ref. ll|, a deviation from 
the ideal scaling between adjacent energy transitions was 
reported, and interpreted as a consequence of electron- 
electron interactions. In contrast, in our data the scaling 
of the transition energies is well preserved in the low en- 
ergy part of the Dirac cone [see Fi^QJb)], while at higher 
energies, we observe deviations due to the non-linearity 
of the dispersion relation, as discussed below. 

A simple theoretical model has been developed to 
investigate deviations from the relativistic Dirac case 
within the tight-binding model with a nearest-neighbor 
(n.n.) hopping term t ~ 3 eV on a honeycomb lattice [l7j ]. 
Next-nearest-neighbor (n.n.n.) hopping t', between sites 
on the same sublattice, is also included, with t' /t ~ 0.1 
[l8j |. The model considers a single graphene layer, which 
is a valid assumption in the case of almost decoupled 
layers. The energy dispersion for this model may be ob- 
tained from a diagonalization of the 2x2 Hamiltonian 
matrix, which reflects the presence of two triangular sub- 
lattices A and B, 



W(q) 



U(q) 



h'(q) 



(1) 



ity of multi-layer graphene, as compared to graphite, is 
due to rotational disorder in the stacking, which reduces 
strongly the interlayer coupling by roughly two orders 
of magnitude [l6| . To cover the full spectral range, two 
different experiments have been performed. The unpo- 
larized far-infrared magneto-transmission of the sample 
at T = 1.9 K has been measured using a Fourier Trans- 
form Spectroscopy (FTS). To explore the higher energy 
range the magneto-transmission up to the visible light 
range has been measured at T — 4.2 K using a tungsten 
halogen lamp. Transmission measurements were circu- 
lar polarization resolved and spectra recorded for both 
polarities of the magnetic field. A representative trans- 
mission spectra is shown in Fig. [Ha) . All spectra show 
a number of absorption lines which can be assigned to 
transitions between L_ m (_„) and L n t m ) LLs, where m, n 
enumerates the Landau levels. In the experiment, we 
observed all transitions from L_i — > L2 (£-2 — > L\) to 
L-13 — * £12 (L_i2 — > £13)- 

The energetic position of the observed absorption lines 
is plotted as a function of the square root of the magnetic 
field in FigQJb). At low energies, the positions of the op- 
tical transitions follow the theoretical (linear) prediction 
for Dirac particles [dashed lines in FigQJb)]. For energies 
above ~ 500 meV a deviation from the predicted linear 
Dirac dispersion starts to be observed (this can be seen 
more clearly in Fig. [2|). Recently Jiang et al. ll| pre- 
sented FIR spectra of a single layer of exfoliated graphene 



with h(q) = — t exp(iq • a 3 ) and h'(q) = 
2t'^2j =1 cos(q • Tj). Here, the vectors ai = a(\/3e x + 
e y )/2, &2 = 5(— ^/ie x + e y )/2, and a3 = — ae y indicate 
the coordinates of n.n. carbon atoms, with a distance 
a = 0.14 nm, and t\ — V3ae x , T2 — V3a(e x +-\Z3e y )/2, 
and T3 = v / 3a(— + \/3e y )/2 those between n.n.n. 

In order to obtain the low-energy spectrum of the dis- 
persion, one expands h(q) and h'(q) around the K and 
K' points at the edges of the first Brillouin zone (BZ), 
characterized by the wave vectors ±K = ±(47r/3v / 3a)e a; . 
An expansion in k = q =p K, up to third order yields 

H+,k) = ^(V^fc* 2 -^|k| 2 fc) ( 2 ) 

M-.k) = -nvp(k*+~^k 2 -^M 2 k*y (3) 

for the K (a — +) and K' (a = — ) points, respectively, 
where we have used the complex notation k = k x + ik y 
and Vp = 3ta/2h. Furthermore, we have introduced the 
phenomenological parameters Wi and W2 , in order to ac- 
count for corrections beyond the simplest tight-binding 
model (l9j |. which yields w± = w^, = 1. For the n.n.n. 
term, expanded to lowest non-trivial order around K and 
K' one obtains 
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and thus the total energy dispersion, taking into account 
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both higher-order band corrections and n.n.n. hopping, 
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where a = + denotes the conduction and a = — the va- 
lence band, and tan 0k = k y /k x . We have subtracted the 
unimportant constant — 3t' , redefining the zero-energy 
position. The cosine terms in Eq. ([5]) indicate that the 
energy dispersion becomes anisotropic (trigonal warping) 
[HI. One clearly notices from Eq. ([5]) that n.n.n. hop- 
ping breaks the particle-hole symmetry but leaves the 
Fermi velocity unchanged. 

In order to account for the magnetic field, we use the 
Peierls substitution, which consists of replacing the wave 
vector k by a momentum operator in the continuum mini- 
mally coupled to the vector potential A, k — * II = p+eA 
[20| . The operator II may be expressed in terms of 
harmonic-oscillator ladder operators, with [a, a'] = 1, 
and the Peierls substitution thus reads 



and 
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Here, l B = yJKJeB = 26/y r B\T\ nm is the magnetic 
length, which is large in comparison with a, and the 
above corrections to the linear energy dispersion are gov- 
erned, in the presence of a magnetic field, by the small 
parameter a/Is- The substitution ([S]), together with Eqs. 
@-(H]), allows one to calculate the energies of the rela- 
tivistic LLs, which, in the absence of the trigonal- warping 
terms, read 
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Here, we have defined 7 = \/2Hvf/Ib and neglected 
terms due to the order of the operators a and a' when 
using the substitution ([6]). This is justified in the large-n 
(semiclassical) limit. In order to account for trigonal- 
warping at leading order, we use perturbation theory, 
which is justified because cl/Ib -C 1. There is no con- 
tribution at first order since (n|a^ 3 |n) = due to the 
orthogonality of the eigenstates (n\n'} — <5 n , n '. At second 
order, one obtains — 7 2 w 2 (a/^s) 2 [3n(n + 1) + 2]/8, which 
needs to be added to the l.h.s. in Eq. |(7J). The fact that 
trigonal warping is manifest only at order {cl/Ib) 2 is due 
to the magnetic field, which averages to zero the cos 30k 
term in Eq. ([5]) when summing over the angle 0k- 

One finally obtains, in the large-n limit, where these 
corrections become relevant, the energies of the relativis- 
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FIG. 2: (color online) The deviation from linearity AE (from 
the data in Fig. [TJb), with the same colors) for the different 
transitions as a function of the energy A°. The solid lines 
is the result of the theoretical calculation as described in the 
text. 



tic LLs for both valleys, K and K', 
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where 0{n ) stands for corrections of order unity, and 
we have defined w 2 = {w 2 + 2w|)/3, which may be mea- 
sured experimentally. The LL structure in the presence 
of n.n.n. hopping has been discussed before in Ref. 21 1. 
We have checked the above result within the semiclassi- 
cal Onsager quantization scheme, and a comparison with 
a numerical solution of the Harper equation on the hon- 
eycomb lattice shows excellent agreement even at small 
values of n [22| ■ 

In our experimental study, we compare the deviation 
AE^ 1 = A° — A„ between the interband-transition en- 
ergies = j{y/n + 1 + yfn) of the ideal case of Dirac 
electrons with linear dispersion and the measured transi- 
tions A+ (for — n — > (n+1)) and A^ (for —(n + 1) — > n), 
as a function of A° . In Fig. [21 we compare these devia- 
tions to the theoretical expectations 
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obtained from Eq. ([5]) in the large-n limit and A^ = 
±(e±,n+i — £=f,«)- Note that A° cx \f~Bn when n> 1. 
The solid lines in Fig. [2] show the theoretical result ([9]) 
with a fitting parameter w = 2.8 (compared to w = 1 in 
the simplest tight binding model). This somewhat large 
value indicates that although the tight-binding model 
yields the correct functional form and order of magni- 
tude of band corrections, it underestimates the strength 
of these corrections, in particular the effect of the trigonal 
warping. In order to account for this enhanced value in a 
theoretical model, one would need to include corrections 
due to the overlap of the atomic wave-functions on the 
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FIG. 3: (color online) (a) The polarization selection rules for 
optical transitions in graphene. (b) The differential transmis- 
sion spectra for both circular polarization of the light. 



different lattice sites, larger distance hopping, and pos- 
sibly interaction effects. One may also speculate that, 
although the graphene layers are only weakly coupled, 
the remaining interlayer coupling might play a role [23| . 
Indeed, trigonal warping in bilayer graphene and graphite 
is dominated by interlayer hopping, which could be on 
the same order of magnitude as the abovementioned dis- 
persion corrections [ij, [25[ . 

Eq. © shows that the LLs are not electron-hole sym- 
metric due to n.n.n. hopping. The resulting asym- 
metry A = |A~ — A+| in optical LL transitions is 
A = 3V2-yt'a/tl B ~ 0. 08 B[T] meV, which is independent 
of n. At fields as high as 30 T, one therefore expects an 
electron-hole asymmetry on the order of 2.5 meV, which 
roughly corresponds to the thickness of the theoretical 
curve in Fig. [51 The effect is thus beyond the resolution 
of our experimental data. 

Additional confirmation of the small electron-hole 
asymmetry in the LL transitions can be seen using a 
polarization- resolved experiment. The polarization se- 
lection rules for optical transition in graphene are shown 
schematically in Fig. [3fa). In Fig.^b) we present trans- 
mission spectra measured for both polarizations at differ- 
ent magnetic field values. No significant differences be- 
tween the positions of the absorption lines can be seen for 
the different polarization suggesting that there is no ob- 
servable asymmetry between the electron and hole cones. 

Symmetry breaking in a gated single sheet of graphene 
has been reported recently by Deacon et al. [13] • In their 
cyclotron resonance measurements, the asymmetry is at- 
tributed to n.n. wave-function overlap corrections char- 
acterized by the overlap integral sq [26] . In this case 
the LL transition asymmetry is A' = ?>\/2 r y sqcl / 1 b oc B, 
which shows that so plays a role similar to t'/t, even 
though the two types of asymmetry have different micro- 
scopic origins. Deacon et al. estimated the strength of 
the asymmetry to be A' ~ 5 meV at B ~ 9 T from the 
— ► 1 and —1^0 transitions [T3], which is in between 
5 and 7 times larger than theoretical estimates, depend- 
ing on whether one takes sq = 0.129 [26[ or t'/t ~ 0.1 



[18| with t ss 3 eV. For a field of 32 T, this would yield 
an asymmetry of the order of 18 meV, which should be 
visible, but is clearly not observed in the 32T spectra in 
Fig.EIb). 

In conclusion, we have probed the high-energy range 
(< 1.25 eV) of the Dirac cone in multi-layer graphene 
and observed a significant deviation from the linear dis- 
persion for massless Dirac fermions. A theoretical model 
which includes higher order band terms gives good agree- 
ment with experiment. No electron-hole asymmetry in 
interband-LL excitations is observed, in agreement with 
our theoretical description, where this asymmetry plays a 
minor role as compared to trigonal warping of the Fermi 
surface and higher order band corrections. 
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